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ABSTRACT 

Context. The nearby K2 V star e Eridani hosts one known inner planet, an outer Kuiper belt analog, and an inner 
disk of warm dust. Spitzer/IRS measurements indicate that the warm dust is present at distances as close as a few AU 
from the star. Its origin is puzzling, since an "asteroid belt" that could produce this dust would be unstable because of 
the known inner planet. 

Aims. Here we test a hypothesis that the observed warm dust is generated by collisions in the outer belt and is trans- 
ported inward by Poynting-Robertson drag and strong stellar winds. 

Methods. We simulated a steady-state distribution of dust particles outside 10 AU with a collisional code and in the 
inner region (r < 10 AU) with single-particle numerical integrations. By assuming homogeneous spherical dust grains 
composed of ice and astrosilicate, we calculated the thermal emission of the dust and compared it with observations. 
We investigated two different orbital configurations for the inner planet inferred from radial velocity measurements, 
one with a highly eccentric orbit of e = 0.7 and another one with a moderate eccentricity of e = 0.25. We also produced 
a simulation without a planet. 

Results. Our models can reproduce the shape and magnitude of the observed spectral energy distribution from mid- 
infrared to submillimeter wavelengths, as well as the Spitzer/MIPS radial brightness profiles. The best-fit dust compo- 
sition includes both water ice and silicates. The results are similar for the two possible planetary orbits and without a 
planet. 

Conclusions. The observed warm dust in the s Eridani system can indeed stem from the outer "Kuiper belt" and be 
transported inward by Poynting-Robertson and stellar wind drag. The inner planet has little effect on the distribution of 
dust, so that the planetary orbit could not be constrained. Reasonable agreement between the model and observations 
can only be achieved by relaxing the assumption of purely silicate dust and assuming a mixture of silicate and water 
ice in comparable amounts. 



Key words, planetary systems: formation 
- zodiacal dust 



1. Introduction 



circumstellar matter - stars: individual: e Eridani - planet-disk interactions 



The nearby (~ 3.2 pc) K2 V star e Eridani (HD 22049, 
HIP 16537 HR 1084), with an age of & 1 Gyr 
dSaffe et al.1 [20051 iDi Folco et all [2001 iSong et all [20001 
ISoderblom fc Dappen 1989T) . has a ring of cold dust 
at ~ 65 AU seen i n re solved sub-millimeter images 
([Greaves et al.l 119981 |2005[) . which encompasses an m- 
ner disk of warm dust revealed by Spitzer/MIPS 
(jBackman et al.l I2009D. The s t ar is orbited by a radial 
velocity planet ([Hatzes et al.l l2000t ) with a semimajor 
axis of 3.4 AU. Another outer planet may orbit near ~ 
40 AU, producing the inner cavity and clumpy structure 
in the outer ring (iLiou fc Zooklll999t lOzernov et al.| [2000t 



Here, we check the possibility that the source of the 
warm dust is the outer ring from which dust grains could 
be transported inward by Poynting-Robertson drag and 
stellar wind. The importance of the latter for debris disks 
around late-type, low -mass stars was first pointed out by 
iPlavchan et al.l ()2005[ ). and it is well know n that e Eridani 
does have strong winds ([Wood et al.ll2002ri . 



iQuillen fc Thorndikdl2002HDeller fc Maddisonl 1200511 . The 
excess emission at A > 15 pm in a Spitzer/IRS-spectrum 
(jBackman et al.l I2009D indicates that there is warm dust 
close to the star, at a few AU (inset in Fig. Q]). Its origin 
is unknown, as an inner "asteroid belt" that could produce 
this dust would be dynamically unst able because of the 
known inner planet ([Brogi et al.ll2009l) . 
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It is convenient to divide the entire system into three 
regions: the outer "Kuiper belt" (55-90 AU), the intermedi- 
ate zone (10-55 AU), and the inner region (inside 10 AU), 
see Fig. [T] In Sect. [2] we describe the modeling setup. In 
Sect. [3] we model the dust production in the outer ring and 
its transport through the intermediate region. Section0]de- 
scribes simulations of the inner system. Section [5] presents 
the spectral energy distribution (SED) modeling and pro- 
vides an additional check for connection between the inner 
system and the outer parent ring (Sect. [53]). Section [6] fo- 
cuses on the surface brightness profiles. Section [7] contains 
our conclusions. 
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Fig. 1. A schematic view of the e Eridani system's archi- 
tecture. The outer ring is the region where the dust is pro- 
duced by parent planetesimals; the intermediate zone is the 
one where it is transported inward by drag forces, possibly 
interacting with a presumed outer planet; transport contin- 
ues through the inner region where dust interacts with the 
known inner planet. Two possible orbits of the inner planet 
are shown. The outer part of the sketch (> 10 AU) is not to 
scale. Inset: The observed SED. The IRS spectrum (dots) 
stems from dust in the inner region and exhibits a charac- 
teristic "plateau" ( "shoulder" ) at Aw 20-30 pm. The main 
part of the SED with a maximum at A k 70-80 pm, well 
probed by several photometry points (symbols with error 
bars), derives from the outer and intermediate regions. 



2. Model setup 

2.1. Method 

Our model includes gravitational forces from the star and 
an inner planet, radiation and stellar wind pressure, drag 
forces induced by both s tellar photons and stellar wind par- 
ticles (jBurns et al.lll979t ). as well as collisions. Dust produc- 
tion in the outer region and dust transport through the in- 
termediate region are modeled with a statistical collisional 
code. To study the dust in the inner region, we need to han- 
dle dust interactions with the inner planet. These cannot 
be treated by the collisional code, so we model the inner 
region by collisionlcss numerical integration. 



2.2. Stellar properties 

We assumed a ste llar mass of M.+ = 0.83A^© 
Benedict et al.l 20061) and a luminosity of L+ = O.32L 



..Pi Folco et al. I l2007t). ? or the stellar spectrum we used a 
NextGen model (jHauschildt et al.lll999n with an effective 
temperature of 5200 K, logq = 4.5, solar metalli city, and 
stellar radius R+ = 0.735i? Q (|Di Folco et al.ll2007l) . 



70% ice + 30% astrosil 
50% ice + 50% astrosil 
100% astrosil 
blow-out limit 




10" 

grain size [um] 

Fig. 2. The /3 ratio of two silicate - water ice mixtures, 
compared to a pure silicate composition, as a function of 
size for e Eri. The bars show the size bins used in the inner 
disk simulations (Sect.|5|). The dashed horizontal line shows 
the dynamical blowout limit of (3 = 0.5. 



2.3. Dust grain properties 

The knee in the IRS spectrum at ~ 20 pm (inset in Fig. [T]) 
is reminiscent of a classical silicate feature. Since the ex- 
act composition of those sili cates is not known, w e have 
chosen astronomical silicate (jLaor &: Draind I1993I ) (pd = 
3.5gcm~ 3 ). On the other hand, by analogy with the sur- 
face composition of Kuip er belt objects in the solar system 
(e.g.. lBarucci et al]|2008l ). we may expect many additional 
species such as ices and organic solids. In particular, it is 
natural to expect water ice to be present, especially given 
that the source of dust is a "Kuiper belt" located very far 
from the star (~ 55-90 AU), and the star itself has a late 
spectral class. Accordingly, we also tried homogeneous mix- 
tures of astr osilicate with 50% an d 70% volume fraction 
of water ice (|Li fc Greenberdll998f) (p d = 1.2gcm- 3 ). The 
bulk density of these ice-silicate mixtures is p — 2.35 gem -3 
and p = 1.89 gcm~ 3 , respectively. The optical constants of 
the mixtures were calculated by effective medium theory 
with the Bruggeman mixing rule. In all three cases (pure 
astrosilicate, two ice-silicate mixtures) the dust grains were 
assumed to be compact spheres. 



2.4. Radiation pressure 

Using the optical constants and adopting Mie theory, 
the radiation pressure efficiency Q Tp , averaged over the 
stellar spectrum, was calculated as a function of size s 
(jBurns et al.lll979HGustafsonlll994h . We then computed the 
radiation pressure to gravity ratio, (Fig. [2]). The resulting 
P(s) was utilized to compute the direct radiation pressure 
and Poynting-Robertson forces. 



2.5. Stellar wind 

The stellar wind was included by a factor /3 sw //3, which is 
the ratio of stellar wind pressure to radiation pressure: 



M 



F„ 



Qrp 



(i) 



where Q sw is the efficiency factor f or stellar wind pressure 
([Burns et al.lll979l: lGustafsonlll994h . We adopted Q sw = 1. 
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Assuming a stellar wind velocity equal to the average so- 
lar wind velocity, v sw ~ 400 km s -1 , and using a mass- 
loss rate of M+ 3OA4 (|Wood et all l2002h . we get 
Asw = 3.69 x 10~ 2 /3. In this estimate (but not in the simula- 
tions) , we set Qrp to unity. Thus direct stellar wind pressure 
is approximately 27 times weaker than radiation pressure. 
However, the stellar wind drag is (/3 SW / (3)(c/v sw ) « 28 times 
stronger than radiation (Poynting-Roberston) drag. 

2.6. Sublimation 

Using an ice-silicate mixture for dust raises a ques- 
tion whether and where the icy portion of the dust 
grains will suffer sublima tion. The ice sublimation tem- 
perature of ~ 100 K ( e.g. iMoro-Martm fe Malhotral l2002t 



Kobavashi et al 



- l2008h is reached at « 10AU (Fig. 
Mukai fc Fechtid (1 1 9831 ) proposed that fluffy and nearly ho- 
mogeneous grains of ice and silicates can produce a core of 
silicate after sublimation of ice. Assuming a grain with a 
volatile icy mantle , instead of a homogeneous sphere of ice 
and silicates, would lead to the same result. In both cases, 
the silicate cores left after ice sublimation continue to drift 
further inward. Within the ice sublimation distance, the 
optical depth w ould decrease by the vo lume fraction of the 
refractories (see lKobavashi et al.ll2008l . their Eq. (35)). The 
shape of the size distribution is not affected by sublimation, 
as long as the volume fraction of ice is independent of the 
grain size. 
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Fig. 3. Equilibrium grain temperature of an 70% ice and 
30% astronomical silicate composition depending on grain 
size s and distance to the star r. The results for the 50% ice 
- 50% silicate mixture are similar. The wavy shape of the 
isotherms is caused by the dependence of t he absorption 
efficie ncy on the grain size (as explained in iKrivov et al.l 
I20081 see their Figs. 2 and 6). 



Thus we expect the entire dust disk to consist of an 
inner silicate disk (r ^ 10 AU) and an outer disk of an ice- 
silicate composition. Accordingly, in integrations described 
in Sect. 0] we assumed dust inside 10 AU to consist of pure 



astrosilicate. Outside 10 AU, we tried all three dust compo- 
sitions described above. 

We fin ally note that a r i ng due to s ublimation as de- 
scribed in iKobavashi et al.l (|2008l I2009T) is not expected, 
because such a ring can only be produced by particles with 
high /3 ratios (/3 ^ 0.1-0.3). Such (3 ratios are barely reached 
in the e Eri system even for pure silicate dust. 

3. Dust in the outer and intermediate region 

3.1. Model 

We used our statistical colli sional code ACE (Analysis 
of Collisional E volution) ijKrivov et al.l 12005k 120061: 
lL5hne et alJ l2008j ) to model the collisional disk be- 
yond 10 AU. This i ncludes the parent ring near 65 AU 
([Greaves et al.l l2005h and the intermediate region inside 
it. The ACE simulation provides us with a rotationally- 
symmetric steady-state dust size distribution at 10 AU, 
which we use later as input for our non-collisional mod- 
els of the inner disk (see Sect. The code is not able 
to treat planetary perturbations. Thus we neglect the pre- 
sumed outer planet, but will qualitatively discuss its possi- 
ble influence in Sects. I3~2l andl5l 

We made three ACE runs, assuming the dust composi- 
tion to be pure astrosilicate, a mixture of ice and astrosil 
in equal parts, and a mixture of 70% ice and 30% astrosil, 
as described above. In handling the collisions, the critical 
specific energy for fragmentation and dispersal QJj is cal- 
culated by the sum of two power laws, 



= A s 



s 
1 m 



3b, 



A B 



3b x 



1 km 



) 



(2) 



where the first and the second terms account for the 
strength and the gravity regimes, respectively. We took the 
values A s = 10 7 erg g~ 1 , A s = 10 7 ergg- 1 3b s = -0.37, 
and 3fe g = 1.38 (cf. iBenz fc Asphaud Il999h . In all ACE 
runs, we assumed that the outer ring has a uniform surface 
density between 55 and 90 AU and that the eccentricities of 
parent planetcsimals range between 0.0 and 0.05, and arbi- 
trarily set the semi-opening angle of the disk to 3°. Next, 
we assumed that the initial size distribution index of solids 
fro m km-sized planet esimals down to dust is q ~ 3.6 (see, 
e.g. lLohne et alj2008t for a justification of this choice). The 
current dust mass was set to M(s < 1mm) w 10 _3 .M®, 
based on previous estima tes from sub-millimeter images 
dGreaves et all 11991 l2005h . 

3.2. Results 

Figure 2] (left) shows the resulting size distribution in an 
evolved disk at several distances from the star. It shows the 
case of a mixture of 70% ice and 30% astrosil, but the results 
for the two other compositions are qualitatively similar. We 
start a discussion of it with the parent ring (55-90 AU) . For 
grains larger than s c ~ 10 pm, the differential size distri- 
bution has a slope q w 3.7. Since Fig. 2] (left) plots the 
optical depth per size decade, this corresponds to a slope 
of g — 3 = —0.7. This is close to what is expected the - 
oretically for collision-dominated disks (|Dohnanvil I1969I ). 
However, from s c « 10 pm down to the smallest grains, the 
size distribution flattens, because the inward drift is faster 
for smaller grains - or, more exactly, for gra i ns with higher 
/3-ratios (see, e.g., IStrubbe fc Cthanel 120061 : IVitense et al.l 
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Fig. 4. Left: Size distribution at different distances in a disk composed of 70% ice and 30% silicate, simulated with 
ACE. Curves depict the normal geometrical optical depth per size decade, as a function of size. The straight lines show 
power laws with q = 3.0 and 3.7 for comparison. (Their height is arbitrary.) Right: Radial profiles of the optical depth. 
Symbols: partial contributions to the optical depth by several selected sizes. Each contribution is normalized to the unit 
size decade around that size. Solid line: the total optical depth of the disk. 



[2010l for a more detailed discussion of this phenomenon). 
The slope calculated analytically for a transport-dominated 
disk with j3 oc s" 1 is q « 2.5. The actual size distribution 
is wavy, and is rather close to q = 3.0, which would corre- 
spond to a uniform distribution of the optical depth across 
the sizes. One reason for that is a nonlinear dependence of 
f3 on the reciprocal of particle size (Fig. [2]). 

Closer to the star, the slope for big grains progressively 
steepens because of their preferential collisional elimina- 
tion. At the same time, the break s c in the size distribution 
moves to smaller sizes. Already at 20 AU, particles larger 
than ~ 10 pm arc almost absent. At 10 AU, the cutoff shifts 
to s ~ 4 um. We argue that a progressive depletion of larger 
grains with decreasing distance will be strengthened fur- 
ther by the alleged outer planet at 40 AU. That planet, 
which presumably sculpts the outer ring, would stop bigger 
grains more effectively than smaller ones, by trapping some 
of them in me an-motion resonances (MMRs) and scatter- 
ing the others ([Liou et aD 1996 : Moro-Martm fc Malhotral 



I2003t iMoro-Martfn fc Malhotral l2005t) . As a result, we do 
not expect any grains larger than about 1-10 um (/3 <, 0.1- 
0.01) throughout most of the intermediate zone and in the 
entire inner region. At the same time, the outer planet is 
not expected to be an obstacle for smaller particles that, as 
discussed below, will be the most important for the mid-IR 
part of the SED. 

From the same ACE run for the outer and intermedi- 
ate regions, we found that r is nearly constant from 55 AU 
down to 10 AU (Fig. [4] right). This optical depth is domi- 
nated by submicron-sized and micron-sized grains. A con- 
stant optical depth inward from the sources is known to 
be a characterist i c feature of transport-dominated disks 
(e.g iBriggsl Il962fc IWvattl [2005h . However, disks with the 
optical depth of r ^ 10~ 5 ai( ? usually thought to be 
collision-dominated (jWvattl [20051 ). and one might not ex- 
pect that even the amplification of the Poynting- Robertson 
drag by stellar wind drag would increase the critical op- 
tical depth (that separates the transport- and collision- 
dominated regimes) by more than an order of magni- 
tude. This raises a natural question: how can it be that 
a disk with the optical depth of ~ 2 x 10 -4 is transport- 
dominated? To find an answer, we note that s c discussed 
above is a critical grain size, for which the collisional life- 



time is equal to the characteristic inward drift timescale. 
Therefore, the disk is collision-dominated at s > s c and 
transport-dominated at s < s c . In the "usual" disks, in- 
creasing r would sooner or later force s c to reach the radi- 
ation pressure blowout limit Sblow At that value of r, the 
entire disk becomes collision-dominated, since there is very 
little material in the disk of size s < Sblow But not in the 
e Eri disk! Here, the blowout limit does not exist (Fig. [2]). 
Therefore, for any reasonable size distribution (q > 3), the 
optical depth is dominated by small particles with s < s c , 
and these are in the transport-dominated regime. That is 
why our simulation shows that the outer and intermediate 
regions of the e Eri disk is dominated by grains roughly in 
the 0.05-1 pm size range and, at those sizes, is nearly colli- 
sionless, despite r of ~ 10 -4 . The same conclusion holds, of 
course, for the inner disk inside 10 AU. For this reason, we 
believe our model of the inner disk in the following section, 
which is obtained via collisionlcss numerical integration, is 
appropriate. 



4. Dust in the inner region 

4.1. Model 

To investigate the behavior of dust in the inner region 
(r < 10 AU) of the e Eri debris disk and its interaction 
with the inner planet, we performed numerical integrations 
of grain trajector ies. We used a Burlisch-Stocr algorithm 
(jPress et al. Il992f ) . Collisions are not considered in our in- 
ner disk model, because, as shown above, they play a minor 
role in the inner disk. 

The dust grains were treated as massless particles, de- 
scribed only by their (5 ratio (or sizes) and their orbital 
elements: semimajor axis a, eccentricity e, inclination i, lon- 
gitude of ascending node f2, argument of pericenter w, and 
mean anomaly M . 

In our simulations we examined the following grain sizes 
s t (i = 1,...5): 0.07, 0.11,0.43, 0.84, and 2.42 pm. These cor- 
respond to the /3 ratios of (3i = 0.1, 0.2, 0.2, 0.1, and 0.03, re- 
spectively (Fig. [2]). Each Si represents a size interval [§i, Si]. 
The limits of these intervals are set to the middle of Si 
and the adjacent sizes Si-i and s^+i (in logarithmic scale). 
Since, due to the low luminosity of the central star, the 
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blowout grain size does not exist, the lower cut-off for the 
particle sizes was arbitrarily set to s\ = 0.05 um. Smaller 
particles are not expected to contribute significantly to the 
SED at mid-IR and longer wavelengths. In addition, we ex- 
pect that various erosive effects (e.g., plasma sputtering) 
and dynamical effects (e.g., the Lorentz force), which are 
not included in our model, would swiftly eliminate the tini- 
est grains from the system. Particles larger than £5 = 4 um 
are not considered in these simulations, because they are 
absent in this region (cf. Fig. 01 right). 

Since grains with the same /3 ratios experience the same 
force, it was enough to run three simulations (/% = 0.2, 
0.1, and 0.03) to cover all five grain sizes. We used 10 000 
particles for each value of f3. The grains were placed in or- 
bits with initial eccentricities from to 0.3 and an initial 
scmimajor axis of 20 Alfl, assuming that they have passed 
the expected outer planet and are out of the range of its 
perturbation. The rather high initial eccentricities of up 
to 0.3 were taken, because these are expected to be in- 
creased by perturbations of the outer planet when passing 
through its orbit. The initial inclinations were uniformly 
distributed between 0°and 10°, respectively. The angular 
elements (£l,u>,M) were all uniformly distributed between 
0°and 360°. 

For the inner planet we used two different orbital con- 
figurations. Th e se con figur ations fit the observ ations from 
IBenedict et"aH (|2006l ) and iButler et all (|2006t) . in the fol- 
lowing called cases 'A' and 'B', respectively. In both cases 
we set the planet's semimajor axis to 3.4 AU and as- 
sumed its orbit to be co-planar with the disk's mid-plane. 
IBenedict et al.l (|2006l ) estimate the planet's mass and ec- 
centricity to M. p d = 0.78A^j up and ep A = 0.7, while 
IButler et all (|2006f ) give M Pb = 1.06X Jup and ep B =0.25. 
We also produced a simulation without an inner planet. A 
summary of the initial parameters for the simulations is 
given in Table [TJ 

4.2. Results 

From our simulations we created a radial profile of the nor- 
mal geometrical optical depth, assuming a steady-state dis- 
tribution (Fig. [5]). It illustrates that the inner planet inter- 
cepts a fraction of dust on its way inward by capturing the 
grains in MMRs and/or scattering those grains that pass 
the planetary orbit. This results in an inner gap around and 
inside the planet orbit. The gap is deeper for grains with 
smaller (3 ratios, because they drift inward more slowly, 
which increases the probability of a resonance trapping or 
a close encounter with the planet. Another conclusion from 
Fig. [5] is that the radial profiles in the cases 'A' and 'B' 
do not differ much, although the eccentricities of planetary 
orbits are quite different. 

5. Spectral energy distribution 

5.1. SED from the outer and intermediate regions 

Having calculated the distribution of dust particles with dif- 
ferent j3 ratios in the s Eridani system, we need to calculate 



Table 1. Initial parameters of the inner disk simulations 



star 


M* 

L* 

M* 


O.83A4 
0.32L Q 

30Mq = 6 x lO^'Meyr" 1 
400 kms -1 


planet 


A B 


Mp 
a P 

Cp 

ip 

Mp , fip , ujp 


1.78 x 10~ 3 .M* 2.44 x 10~ 3 .M* 
3.4 AU 3.4 AU 
0.7 0.25 
0.0° 0.0° 
0.0° 0.0° 


dust 


P 

do 

eo 
io 

Mo,fi ,wo 


0.2, 0.1, 0.03, 0.01 
20 AU 
[0.0, 0.3] 
[0°, 10°] 
[0°,360°] 



1 We use 20 AU instead of 10 AU because, if we placed the 
particles at 10 AU with eccentricities from to 0.3, the initial 
distances would be distributed between 7 and 13 AU, causing an 
unwanted "boundary effect" of decreased optical depth in that 
distance range. 



A range of [a, b] indicates that the parameters are equally 
distributed between a and b, and randomly chosen. 



the thermal emission of the dust for comparison with the 
observed SED. We start with the SED of the outer and in- 
termediate region together. We have calculated it, using the 
output of the ACE runs described in Sect. [3] Specifically, 
coupled radial and size distributions of dust between 10 and 
90 AU shown in Fig. 0] were used. 

The results are shown in Fig. R)! Overplotted are the data 
points, all taken from Fig. 7 of iBackman et all (|2009l ): see 
also their Table 1. Clearly, the model with 100% astrosili- 
cate contents is too "warm" . Notably, the SED starts to rise 
towards the main maximum too early (at « 25 um instead 
of ps 35 jim) , which is inconsistent with the "plateau" of the 
IRS spectrum. Including 50% ice improves the model, but 
the SED still rises too early. The best match of the data 
points is achieved with a 30% astrosil and 70% ice mixture, 
which will therefore be used in the rest of the paper. Still, 
even for that mixture, the main part of the SED (i.e. the 
combined contribution of the outer ring and the intermedi- 
ate region) predicted by the collisional model is somewhat 
"warmer" than it should be. We discuss this slight discrep- 
ancy in Sect. 15.31 



5.2. SED from the inner region 

We now consider an SED coming from dust in the inner 
disk. The optical constants, density of the dust, and the 
stellar spectrum described in Sect. 12.31 have been used to 
calculate the "partial" SEDs from each grain size bin. To 
sum up the contributions from different bins, we assumed 
a power law A4(s)ds oc s 3 ~ q ds. There are two parame- 
ters that allow us to adjust the resulting SED of the inner 
region to the IRS spectrum. One is the slope q that pro- 
vides a "weighting" of relative contributions of different- 
sized grains into the resulting SED. In accord with our re- 
sults in Sect. 13.21 we take q = 3 and only consider grains 
with s Ss 2 pm (four lowest size bins). Another parameter 
should characterize the absolute amount of dust in the in- 



5 



Reidemeister et al.: e Eridani 



Planet A: a=3.4AU, e=0.70 



10"' 



10 c 



0.05-0.09nm 
0.09-0.22|im 
0.22-0.60jim 
0.60-1 .42nm 




4 6 
distance [AU] 

Planet B: a=3.4AU, e=0.25 



10 



10"' 



0.05-0.09jim 
0.09-0.22nm 
0.22-0.60|im 
0.60-1.42nm 




4 6 
distance [AU] 

without planet 



10"' 



« 10- 



10"' 



0.05-0.09nm 
0.09-0.22|im 
0.22-0.60|im 
0.60-1.42nm 




4 6 
distance [AU] 



Fig. 5. Radial profile of the normal geometrical optical 
depth produced by different grain size bins for the planet 
'A' (top), 'B' (middle), and without a planet (bottom). A 
radial bin width is Ar = 0.25 AU. Vertical dashed lines 
show periastron and apastron of the planet orbit. A size 
distribution with a slope of q = 3.0 and the dust mass of 
M dust w 8.5 x 10" 8 A4 e (case 'A'), 9.9 x 10" 8 X e (case 
'B'), and 5.6 x 10 _8 A / I© (without a planet) were used for 
the vertical scaling of the curves. 



ner region. (This could be, for instance, the total dust mass 
in the inner region.) This parameter would determine the 
overall height of the SED. We varied the dust mass until 
the resulting model SED fits the IRS spectrum best. The 
best results were achieved with the total dust mass in the 
inner region of A^dust ~ 8.5 x 10^ 8 A4® in the case 'A', 
9.9 x 10~ 8 A^e in the case 'B', and 5.6 x 10~ 8 A^© without 
a planet. The absolute height of optical depth profiles from 
various size bins shown in Fig. [5] corresponds to the same 
dust masses and the same slope q = 3.0. 
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Fig. 6. The SED from dust outside 10 AU in the e Eri dust 
disk. Symbols with error bars are data points and scattered 
dots are the IRS spectrum. Lines are model SEDs for dif- 
ferent dust compositions: 100% astrosil (dotted), 50% ice 
and 50% astrosil (dashed), and 70% ice and 30% astrosil 
(solid). 

Figure [7] shows the contribution of the different grain 
sizes to the SED and their total, we now use a logarithmic 
scale. Although the two planet orbits are quite different, the 
influence of the planet on the SED is rather minor, because 
the radial distributions of dust are similar in both cases. 

A consistency check that we made was to compare 
the model predictions with the results of intcrferometric 
measurements with C HARA array in the X-band (2.2 pm) 
(|Di Folco et alJl2007h . They set the upper limit of the frac- 
tional excess emission of the inner debris disk to 6 x 10 ~ 3 
(3cr upper limit). With the photospheric flux of 120 Jy at 
A = 2.2 pm, this translates to an excess of Ss 720 mJy. This 
value includes both thermal emission and scattered light. 
The integrated surface brightness of the 2.2 pm radial ther- 
mal emission profile, convolved with the CHARA trans- 
mission profile, generates a total excess of just 25.7 mJy, 
14.5 mJy, and 15.4 mJy for the cases 'A', 'B', and without 
an inner planet, respectively. Even if we took scattered light 
into account, which we estimate to contribute ks 3.4 times 
more than the thermal emission at that wavelength, our 
model would be consistent with non-detection of dust with 
CHARA. 

5.3. SED from the entire disk 

We now assemble the SED produced by the entire disk. To 
this end, we summed up the SEDs of the inner region and 
of the region outside 10 AU presented in Figs. [6] and 
respectively. Figure [5] shows the complete SED. It is in 
a reasonable, although not a perfect, agreement with the 
observations. In particular, the maximum of the modeled 
SED, while reproducing the data points within their er- 
ror bars, appears to lie at a slightly shorter wavelength 
than the one suggested by the data points. A likely reason 
for this discrepancy is that our collisional simulation does 
not take into account elimination of particles in the size 
range from ~ 1 to ~ 100 pm by the alleged outer planet, 
as explained above. Excluding these particles from the in- 
termediate region 10-55 AU would reduce emission in the 
35-70 pm wavelength range, shifting the maximum of the 
SED to a longer wavelength. In addition, the main part of 
the SED can be made "colder" by varying diverse param- 
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Fig. 7. The observed Spitzcr/IRS spectrum 

(jBackman et al.l l2009h (small dots) and the modeled 
SED of the inner region. Top: planet 'A', middle: planet 
'B', bottom: without a planet. Thin lines: contributions 
of the different size bins to the SED, thick line: the total 
emission from all sizes. The size distribution slope q and 
the dust mass A^dust are the same as in Fig. [5] 
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Fig. 8. The entire SED of the e Eri dust disk. Symbols with 
error bars are data points, small dots is the IRS spectrum. 
Thick line in each panel is the excess emission predicted 
by the models, and thin lines show contributions from the 
inner region (< 10 AU) and the intermediate+outer disk. 
Three panels are for the planet configurations 'A' (top), 'B' 
(middle), and without a planet (bottom). 



eters of the collisional simulation, many of which are not 
at all or are poorly constrained. These include the eccen- 
tricity distribution of planctcsimals, the opening angle of 
the planctesimal disk, as well as the mechanical strength of 
solids. Such a search for the best fit would, however, be very 
demanding computationally. We deem the fit presented in 
Fig. [8] sufficiently good to demonstrate that our scenario, 
in which inner warm dust is produced in the outer ring, is 
feasible. 



5.4. Connecting outer and inner regions 

To make sure that the amount of dust required to reproduce 
the IRS spectrum is consistent with the amount of dust that 
could be supplied to the inner disk from the outer parent 
belt, we now make an important consistency check. From 
calculations of the dust production in the outer ring and 
its transport through the intermediate region, we know the 
optical depth at 10 AU. Figure 0] (right) suggests that the 
optical depth per size decade, created by grains with sizes 
of 0.07, 0.11, 0.43, and 0.84 urn, is about (1-2) xlO -4 . These 
values have to be compared with optical depths, sufficient 
to reproduce the IRS spectrum. From Fig. [5] outside the 
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planetary region we read out the values (2-5) xl0~ 5 for 
the four lowest size bins (which are centered on the same 
four grain sizes). Taking into account the actual width of 
the four lowest size bins used in the modeling of the inner 
system, which is w 0.4 dex, we would need an optical depth 
per size decade of (5-10) xl0~ 5 . Thus the optical depth 
that is supplied by grains transported from outside matches 
the optical depth required to account for the IRS spectrum 
within a factor of two. 

6. Surface brightness profiles 

The model developed in this paper was designed to explain 
the Spitzer/IRS spectrum of the system. Besides this, we 
have made sure that it reproduces the entire SED probed 
with many instruments at the far-IR and sub-mm wave- 
lengths. We now provide a comparison with other measure- 
ments that we have not considered before. 

The most important information comes from spatially 
resolved images. Spitzer/MIPS observations yielded bright- 
ness profiles in a ll three wavelength ba nds centered on 24, 
70, and 160 pm ijBackman et alJl2009ft . We calculated the 
brightness profiles of thermal emission with our model. 
In doing so, we included both the outer+intermediatc 
disk and the inner one. The resulting brightness profiles 
were convolved wi th the instrumental PSF (see Sect. 2 in 
iMiiller et alj|2010l for the algorithm used) and compared 
them with observed profiles. 

The results are presented in Fig. [9l We only show case 
'A', since the profiles in case 'B' and without a planet 
are very similar. At all three wavelengths, the modeled 
brightness monotonically increases toward the star, as does 
the observed brightness. Both the slopes and the absolute 
brightness level are in good agreement with observations. 
The only exception is the modeled 70 um profile (Fig. [9] 
middle). It is flatter than the observed one, predicting the 
brightness inside 30 AU correctly but overestimating the 
emission in most of the intermediate and the parent ring 
regions by a factor of two. The reasons for this deviation 
are probably the same as those discussed in Sect. 15.31 First, 
if a presumed outer planet at rj 40 AU efficiently eliminates 
particles in a size range from several to several tens of mi- 
crometers (which is not taken into account in our ACE 
simulations) , this will decrease the 70 pm brightness in the 
intcrmedatc region. Second, it should be possible to slightly 
decrease the 70 pm emission in the parent ring region by 
varying poorly known parameters in the collisional simula- 
tion. 

We also checked the radial brightness profile at 
850 p m and compared it with JCMT/SCUBA observa- 
tions (|Greaves et al.| [l998. 2005j). After convolution with a 
Gaus sian PSF of a = 5", the resulting profile is consistent 
with iGreaves et al.l (|2005l their Fig. 2), showing a broad 
ring around w 60 AU and a resolved central cavity. 

7. Conclusions 

In this paper we show that in the nearby system e Eridani 
the Spitzer/IRS excess emission at A ~ 15-30 pm can be 
caused by dust that is produced in the known outer dust 
ring and that streams inward due to interaction with strong 
stellar winds. 

By running a collisional code, we simulated the dust pro- 
duction in the outer ring between 55 AU and 90 AU with 
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Fig. 9. Azimuthally averaged surface brightness bright- 
ness profiles of the e Eri dust disk (case 'A') at three 
Spitzer/MIPS wavelengths: 24pm (top), 70pm (mid- 
dle), and 160pm (bottom). Points with error bar s : star - 
subtracted MIPS data (Fig. 6 of iBackman et al.l I2009T) . 
Solid lines: modeled profiles after convolution with the 
PSFs. Thin solid lines: the instrumental PSF. 



a dust mass of lO -3 ^® and the subsequent transport of 
the dust inward to 10 AU. We then employed single-particle 
numerical integrations to simulate the dust transport fur- 
ther inward through the orbit of the known inner planet. 
The dust in the inner region was found to consist of grains 
smaller than « 2 pm, and the dust mass inside 10 AU was 
estimated to be (6-10) x 10~ 8 .M©. Combining the results 
of the collisional simulations outside 10 AU and numerical 
integrations inside that distance, we calculated the overall 
SED and radial brightness profiles. This SED is in a rea- 
sonable agreement with the available observational data, 
and it correctly reproduces the shape and the height of the 
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Spitzer/IRS spectrum. Likewise, the brightness profiles are 
consistent with the Spitzer/MIPS data. 

Th e best results are obta i ned with an ice-silicate compo- 
sition (|Laor fc Draind[T99l iLi fc Greenber^[l998h of dust 
outside the ice sublimation distance of « 10 AU, and an in- 
ner disk of non- volatile silicate grains inside that distance. 

With the aid of the modeled spectra and brightness 
profiles, it is not possible to distinguish between the dif- 
ferent orbita l solutions for t h e inn er, r adial velocity plane t 
proposed bv lBenedict et all (|2006l) or iButler et"atl (j2006h . 
Although the planetary orbits they inferred are quite differ- 
ent, both setups yield quite similar radial profiles of dust, 
the SEDs, and the brightness profiles. 

Various kinds of new data on the e Eri system are 
expected soon. Data from JCMT/SCUBA2 (J. Greaves, 
pers. comm.), HERSCHEL/PACS, and /SPIRE should 
shed light on the cold dust. This includes the structure of 
the outer "Kuiper belt" and the intermediate region of the 
disk between 10 and 55 AU. Further more, there is hope to 
better probe the inner warm dust directly, with instruments 
such as the Mid InfraRcd Instrument (MIRI) aboard the 
upcoming James Web Space Telescope (JWST). Likewise, 
there is an ongoing effort to find oute r planets in the system 
by direct imaging (lltoh et all 120061: IMarengo et al.l [20061 



I J anson et afll2007l 120081: IMarengo et al.ll2009t) 
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